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ABSTRACT 


Much has been said in the literature with regard to the problem of 
establishing symbol synchronization in binary baseband digital communication 
systems. By comparison, the literature is virtually devoid of information 
relating to the extraction of symbol sync from multilevel baseband data. With 
the recent interest in multilevel amplitude- shift- keying (MASK) and quadrature 
amplitude- shift- keying {QASK) as signaling techniques for multilevel digital 
communications systems, the problem of providing symbol synchronization in 
the receivers of such systems becomes paramount. 

This paper presents a technique for extracting symbol sync from an 
MASK or QASK signal. The scheme is essentially a generalization of the data 
transition tracking loop (DTTL) which has heretofore been used in PSK systems. 
The performance of the loop is analyzed in terms of its mean- squared symbol 
sync jitter and its effects on the data detection process in MASK and QASK 
systems . 
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DATA-DERIVED SYMBOL SYNCHRONIZATION 
OF MASK AND QASK SIGNALS 

A. INTRODUCTION. 

Much has been said in the literature^ ^ ^ ^with regard to the problem of 

establishing symbol synchronization in binary baseband digital communication 
systems. By comparison, the literature is virtually devoid of information 
relating to the extraction of symbol sync from multilevel baseband data. With 
the recent interest ^ ^ 27 ^in multilevel amplitude- shift -keying (MASK) and 

quadrature amplitude- shift- keying (QASK) as signaling techniques for multi- 
level digital communication systems, the problem of providing symbol synchro- 
nization in the receivers of such systems becomes paramount. 

This paper presents a technique for extracting symbol sync from an MASK 
or QASK signal. The scheme is essentially a generalization of the data transi- 
tion tracking loop (DTTL) ^ 1 which has heretofore been used in PSK 

systems. The performance of the loop is analyzed in terms of its mean- squared 
symbol sync jitter and its effects on the data detection process in MASK and 
QASK systems. 

B. DESCRIPTION OF LOOP OPERATION , 

Consider the symbol synchronizer illustrated in Fig. 1 which is a 
generalization of the DTTL to the case of an MASK input signal. The basic 
pj^ij^ciple of operation of this loop is virtually no different from its binary 
counterpart i n that the upper (in-phase) branch monitors the tran- 

sitions in the input data while the lower (mid-phase) branch produces a measure 
of the lack of synchronization between the input signal and the locally derived 
timing source. The essential differences between the two loops are: (a) the 

hard limiter of the DTTL is replaced by a multilevel quantizer to accommodate 
the multilevel signal input and (b) an additional feed-forward path is provided 
to cancel out the bias voltage introduced by the mid-phase branch even in the 
absence of any synchronization error . Such a bias is not produced in the 
binary DTTL since, in the absence of noise and synchronization error, there 
are only two possible signal levels across a data transition and these levels 
are equal and opposite. 
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The input signal s(t,0 is a random pulse train representing the data 
symbols, i.e. 


s(t,0 = ^ ^ a k p[t - kT - e ] (1) 

k 

where p(t) is the unit rectangular pulse defined by p(t) = 1, 0 < t<T, p(t) = 0; 
all other T, T is the symbol period, a^, the symbol amplitude, takes on 
equally likely values i6 with i = ±1,±3, ■** ±(K-1), and e is the random epoch 
(assumed to be uniformly distributed) to be continuously estimated. The input 
additive noise process, n(t), is assumed to be white Gaussian with a two-sided 
spectral density, N /2 . The sum, y(t), of signal plus noise is input to in-phase 
and midphase filters, so-called because they operate, respectively, across a 
given symbol and between a pair of adjacent symbols. These filters are, in 
principle, matched filters to the input (i.e. integrate and dump circuits), with 
the interval of integration being determined by the estimate e of epoch produced 
by the symbol sync loop. The output y^(t) of the in-phase filter is sampled at 
intervals of T sec. These sample values jy^j are input to a decision device 
and a decision is made corresponding to each multilevel input data symbol. 

The decision device is a multilevel quantizer whose characteristic is given by 

! i&, 

(K- 1)6 , 

-(K- 1 ) 6, 

A correct decision . 



is made whenever ^ ^ corresponding to the transmitted symbol in 

that particular symbol interval. The transition detector then examines two 
adjacent decision a^ ^ and and records an output whenever there is a 
transition in the data, i.e. 


for (i-1) & < x < (i+l)6; 
for (K-2) 6 < x < co 
for - co < x < -(K-2)6. 


i = ± 1 , ±3, ±5, • - - , ±(K-3) 


(2) 
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( 4 ) 


A A 

_ k-1 ~ a k 

A k 2 

We note from (2) and (3) that 1^ ranges over the set of allowable values 
0, ±6, ±26, **• , ±(K-1)6. Thus, the output of the transition detector is both 
a measure of the presence and the magnitude of a data symbol transition. 

The output, y^t), of the mid-phase filter is also sampled at intervals of 
T sec. These sample values jy^^j must be delayed an amount t = (1 - £/2)T 
sec. before multiplication with the appropriate transition detector output I k * 

The quantity |T, 0 < | < 1, is the integration interval of the mid-phase integrate- 
and-dump and is often referred to as its window width . For cases of practical 
interest, where the ratio of two-sided loop bandwidth (to be evaluated later 
on) to data symbol rate <R = l/T is much less than one, the presence of the 
delay t in the loop has negligible effect on its steady- state noise performance. 

We shall denote the output of the delay element by corresponding to the input 

yzk- 

To understand the necessity for the additional feed- forward arm in Fig. 1, 
we examine the sampled output, y^’ o£ the mid- phase filter in the absence of 
additive input noise. If we define the normalized synchronization error by 
\=(€-£)/T, then we can show that 

j k = (l + x ) a k-i + (l - x ) a k = (a k -i + a k>4 - (a k- a k-i )x < 5) 

For the case of binary data, i. e. , = ±6, the presence of a data transition 

results in either = 2X6 or = -2X6 depending upon the direction of the 
transition. Furthermore, if the sync estimation is perfect ( X = 0), then = 0. 
For the case of multilevel data, we note that even for the case of perfect sync, 
a mid-phase filter output voltage still exists, namely J^. = + a^)^/2. 

Thus, in an attempt to remove this bias from the system, we subtract from J^. 
the bias voltage 


[ k = < a k-i + a 


k>t 


(&> 
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The net output of this operation, i.e. - I£. is then multiplied by I to 
produce the error signal e . This error signal is operated upon by a digital 
filter whose output is used to control a timing pulse generator. An estimate, 
e, of the input random epoch is thus formed every T sec. 

In the next section, we consider the steady-state tracking performance 
of the loop in terms of the variance of its timing error (the so-called timing 
jitter). The system parameters which characterize such performance are: 
(1) loop bandwidth- symbol time product W^T, (2) two-sided noise spectral 
density, N Q /2, of the input additive noise, and (3) average signal power S 
where 


S = 6 2 (K 2 - l)/3. 


(?) 


C. SYSTEM ANALYSIS 

We begin with the task of representing our symbol tracker as an equiv- 
alent phase-locked loop whose steady- state performance can be found from a 
well-developed theory (i.e. application of the Fokker-Planck equation). The 
assumptions under which such a representation is valid will be given below. 
Developing the above equivalence relies on finding (a) the average loop S- curve, 
g(\), as a function of the normalized phase error \ = ( e - € )/T, and (b) the 
two-sided spectral density about the origin, S(u>,X), of the equivalent additive 
noise, n^(t), at the output of the loop nonlinearity, g(X). Once having deter- 
mined these quantities, the transition tracking loop may be represented as 
shown in Fig. 2. 

The approach taken in finding g(X) and S(u>, \) follows that of Refs. [7], 

[10]. The assumptions which must be imposed on the system are as follows: 

(A) The input process, €(t), is essentially constant over a large 
number of symbol intervals. 

(B) 6 (t) , the estimate of c(t), is essentially constant over a large 
number of symbol intervals. This fact is a direct consequence 
of a previous assumption that the response of the loop is very 
slow with respect to a symbol interval (i.e. W T « 1). 

Jj 
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It follows from assumptions (A) and (B) that the loop can be analyzed as 
a continuous loop, and that the statistics of the random variable (see Fig. 1) 
can be determined assuming X fixed. In effect, then, what we do is to consider 
many records of the discrete r.v, at fixed X and call the average value of 
this ensemble, g(X), and its spectrum S(w,X). 

Referring to Fig. 1, the output of the in-phase integrate -and- dump 
circuit is 


Yl (kT + i) & y lk = c k + 


( 8 ) 


where 


•kT+e 


c, = 


T 


s(t, e)dt 


'(k- l)T+e 
.kT+€ 


V. = 


k T 


n(t)dt 


/( k-l)T+e 


The output of the mid-phase integrate -and -dump circuit is 


( 9 ) 


y z [( k + i) T + *] = y 2k = b k + |i i 


( 10 ) 


where 


/•(k+g/2)T+J 

b k = I s(t, € )dt 

k-i/2)T+t 

^(k+g/Z).T+e 
=-| n(t)dt 

J( k-^/2)T+c 


(ID 
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Since andv^ are not independent, we wish to express them in terms 
of a new set of random variables: 


\ = N k + M k 


r, = n,' , + m; 

k k- 1 k 


( 12 ) 


where 


N, = — 

1 m 


N k 


•(k-l/2)T+e 


n(t)dt 


(k- 1)T+ C 


*(k-l+|/2)T+€ 


'( k- 1 )T+ ? 


n(t)dt 


M, = b 



kT+£ 


n(t)dt 
r (k-l/2)T+« 


»kT+« 


M' 

k 


n(t)dt 

k'k- 5/2)T + « 


(13) 


with the properties 

N^, are mutually independent for all k, n 

NJ , M' are mutually independent for all k, n 

k n 

N£_, and M^> N are mutually independent for all k, n. 

N, , N 1 and M, , M 1 are mutually independent for all k/n. 

k n k n 

All and their sums are Gaussian random variables with 

zero mean and variances, 



Recalling that y^ 
e k = 






for all k 


- J^(if the mid-phase channel delay is ignored), and 
then combining (1) - (4) and (6) - (14) we get 


(14) 
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e k = w + M k + N k+1 


Q K (c k + M k + N k } + °K (c k+l + M k+1 + N k+l' 


Q K (c k + M k + N k ) - Q K ( c k+ i + M k+ i + N k+ i^ 


(15) 


As previously mentioned, we define the loop nonlinearity, g(\), by 


S (X) = E n, S ( e kl l ) 


where E represents the conditional expectation on X both with respect to the 
in ) s 

noise and signal (symbol sequence). The evaluation of g(X) (for positive X) is 
done in Appendix A, with the result 



E 

m, j, 2 


m(2k) - j(| + X) 


m{2\) + j (A + X) 


erf(j,m,£) - erf(- j, - m,£) 


E 

m, j,£ 


<(m - j) 


-(m + j) 


erf( j , - m, £) - erf( 
0 < X < t/2 


erf(j,m,j£) - erf(-j, - m,£) 


- j , m,l) | 


( 16 ) 


erf(j, - m, £ ) - erf(- j f m,2) 


1 


X < ± 
2 2 


j, 

'Unless otherwise noted, here and in the remainder of the text, sums over m & 
j are for value s m, j=l, 3,5, • • • , (K- 1), while the sum over £ is for values 
j£= 0,2,4, * * * ,(K-2)» It is also easy to show that g(X) is an odd function of X. 
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where ^ is the Neumann number defined by - 1, = 2;k>l, and we have 

introduced the shorthand notation 


erf (j,m, Z) = erf 


^pT(j \ + m( 1 - \) + 



(17) 


with 


R A ST 

R d- 


. R , 

R » £ i 

d 2 


{K - l)/3 


erf < x > 


exp(-t )dt 


(18) 


It is interesting to note that in arriving at { 1 6) , the same result would be 
obtained if the feed-forward signal had been completely eliminated. The con- 
clusion to be drawn from this is that, on the average , the loop S- curve is 
unaffected by the bias previously discussed. Keeping in mind, however, that 
the actual loop operates discretely in that the error signal e^ is updated every 
T sec, then the presence of the feed-forward signal !£. is quite important in 
reducing the instantaneous variations of e^_ about its mean. 

Figures 3 and 4 illustrate the normalized S- curve g^(X) == g(X)/S» defined 
above, with as a parameter, for K = 4 and window widths and £—0.5 
respectively. The limiting S- curve for R d — co is also indicated on these figures. 

Now we examine the problem of finding S(u>,\), the spectrum of the addi- 
tive noise n^(t) = e^t) - g(X), defined by 


S(u>, \) 


= ^{V e k e k +r >> - 


(19) 
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The symbol & denotes the discrete Fourier Transform, i.e. 


CO 

S(w,\) = \ R(m,\) cos wmT (20) 

m=-co 

As for the value of the spectrum as a function of frequency, we note that R(m,\) 
has non-zero value only atm = 0, ±1. Thus, the spectrum consists of the sum 
of a constant and a sinusoidal component with period cJq = 2tt/t. Since as 
before W T T « 1 is assumed, it is sufficient then to consider only the value of 
zero frequency, i.e. S(0,\) and assume a flat spectrum of this value for all to 
of interest. Also, since in a practical system little data degradation due to 
symbol sync error can be tolerated, it is realistic to approximate S(0,K) by its 
value at \=0, viz., S( 0, 0)=R( 0, 0)+2R( 1 , 0) . Such approximations were made 
for the binary data symbol sync loop described in Refs. [7], [10] and were 
shown there to result in negligible performance penalties for system parameters 
of practical interest. The computation of S(0, 0) is a tedious task and is pre- 
sented in Appendix R. The result is 




(ft) 


t R 


2 

K - 1 


K 2 (K 2 - 1) 


exp - £) 2 R^]’ 

E “ E «* [< m ' e) A / 5 i] + — B1 — 1 


K(K - 1) 


2+ |R« 




m l 


( 21 ) 


‘Here the sum over i is for values £-0, ±2, ±4, 


, ±(K- 2) . 
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The normalization in (21) is chosen so that 


limh(O) = 1 (22) 

R cr“ 

Also, for K=2, Eq. (21) reduces to the previously reported result [7] 

£ R d t 

h( 0) = 1 + —— - ~ 

Finally, we note that as in the case of the S- curve calculation, the evaluation of 
the expression for h(0) given in (21) is independent of the presence of the feed- 
forward signal; i.e. the terms involving this signal when averaged over the 
noise and symbol sequence result in zero net contribution. Figure 5, 6 and 7 
illustrate the normalized noise spectrum h( 0) vs. signal-to-noise ratio R^ with 
window width £ as a parameter and K=2,4, 8 respectively. 

D. APPLICATION OF THE FOKKER-PLANCK EQUATION 

The stochastic differential equation which describes the loop illustrated 
in Fig . 2 is 


V R d er£ V r 7 + 


exp(-R d ) 


\]TT 


(23) 


K = - K v F(p) [g(X) + n K (t)] (24) 

where p is the differential operation d/dt and K is the gain of the timing pulse 
generator. In writing the above, we have assumed zero static timing error. 
This is consistent with our previous assumption that e (t) = 0. 

When dealing with equations of the above form, our aim is to find the 

steady-state probability density function (p , d. f, ) of X , p(\) } from which the 

2 

mean-squared timing error, c r^, can be calculated. If we assume that n^(t) i s 
approximately a Gaussian process, then the solution for p(X) can be found via 
the Fokker -Planck method. The validity of this assumption and the manner 
in which the method is applied are discussed in Refs. [7], [1 0] . For sim- 
plicity we will discuss only the first order case (i.e. F(p) = K^) whereupon 
the solution is of the form 


10 
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p(X) = C 1 exp / - 


2R,6,K 
d d g 


g n ( y) d y 


w 4 


where 


W l T 


K K K 
v P 1 
2 


with denoting the two-sided loop bandwidth and K the slope of g (\) at 
X=0; that is (see Appendix C) 


K - y y <- n+K 7 iy 

8 2 < K - *> (n=-K + 3,-iu5, ••• ,-1 K 


( n ^) + 2] er£ k/*I) 

n=l, 3, * • * , K- 1 


n=K+ 1 , K+ 3 , , 2K- 3 




1 ri 

2 V tr 


n=-K+3, -K+5, *** ,-l 


t 1 + iL r i ] exp (- n 2 R d) 


2Z exp f nZR i) 

n=l , 3, • • • , K- 1 


n-K+1 , K+3, **• , 2K - 3 


n -(K - 1) ' 
K 
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Note that 


lim K = 1 (28) 

R d^“ s 

Plotted in Figs. 8, 9, and 10 are K as given by (27) vs. R , with | as a 

§ 

parameter andK=2,4,8 respectively- The normalized mean- squared sync 
error <r 2 is then given by 


2 


1/2 


-l/2 


X p(X)d\ 


(29) 


For large R^& d (the case of practical interest), we can approximate 


g n (X) b Y K Xin Eq. (25). Making this substitution gives 


p(\) = C 2 exp 


2R d 6 d K g 

“pr 




(30) 


Furthermore, since p(\) now has significant values only on a small interval of 
X relative to one -half, then the finite integral in (29) can be approximated by an 
integral from zero to infinity. Thus, by inspection of (30) we have that for 
large R d & d , 


s Z h (°l.y (31 

2R d 6 d K g 

If now in addition R^ is large, then from (22), (28), and (31) we have asymptot- 
ically that 


2 s 5 

2R.6 1 
d d 


(32) 
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which is independent of K and hence is the identical resalt as that obtained for 
the binary symbol sync loop [l], [Yj, [lo]. Plotted in Figs. 11, 12, and 13 is 
the rms symbol sync jitter <r^ [as computed from (31)] vs R^ with 6^ as a param- 
eter, 6=1, and K = 2, 4, 8 respectively. We see from these results that the 
performance of the synchronizer (insofar as this measure is concerned) is virtu- 
ally independent of the number of signal levels K. 

E. ERROR PROBABILITY PERFORMANCE 

This section deals with the effect of symbol sync error on system error 
probability performance. Since the optimum detector for an MASK signal set 
is a matched filter followed by a K-level quantizer, we are interested primarily 
in examining the degradation in error probability performance of this detector 
when the symbol synchronization reference is noisy. The first step in the pre- 
sentation will be to arrive at the error probability of the MASK detector condi - 
tioned on a symbol sync error. This conditional error probability can then be 
averaged over the p.d.f. p(\) of (25) or (30) to yield the average error proba- 
bility performance of the receiver. When compared with a source of perfect 
symbol sync, the degradation in performance due to symbol sync error can be 
assessed. 

Consider the situation where the symbol sync is in error by e = \T sec. 

The matched filter statistic on which a decision is made for the nth transmitted 
symbol now depends on the received data in the time interval (n - 1) T + f < t 
< nT + c where c is the sync estimate obtained from the symbol synchronizer. 
Stated mathematically, 

pT+( j 

| y(t) p[t-(n-l)T-e] dt> (33) 

(n-1 )T + t \ 

where Q K (x) is the K-level quantizer characteristic defined in (7), and 
y(t) = s (t , € ) + n(t) is the matched filter input. Following an approach analogous 
to that given in [l] for binary data, we arrive at an expression for the condi- 
tional error probability P E (\) of an MASK input (see Appendix D), that is 
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p E <"> = ^ 


E 


E 


e rfc 


j=±l, ±3, • ♦ • , ±(K-1) 11 = 0 , ±2, ±4, • • • , ±(K-2) 



R 


i - <« - u (i - |x| ) - j|M 




(34) 


For K = 2, (34) reduces to the well-known result for NRZ data 

P e (M = \ erfc (fc) + 1 erfc (1 - 4 |x|>] (35) 


where 


erfc(x) ^ 1 - erf (x) 


The average error probability is determined from 


1/2 

E ~ / P E ( X ) P(M dx 

- 1/2 


(36) 


(37) 


Assuming large R^6^, then substituting (30) and (34) into (37) and performing 
numerical integration on a digital computer provides the graphical results of 
Figs. 14, 15. In these figures the average error probability is plotted versus 

with 6^ as a parameter for £ - 1 and K = 2, 4 respectively. We see from these 
figures that a 6^ of 250 (for K = 2) and a 5^ of 100 (for K = 4) are sufficient to 
reduce the signal-to-noise ratio degradation (relative to ideal performance) to 
around 0. 1 dB, In a practical system, 6^ would be more of the order of 
thousands. 


F. MODIFICATION OF THE MASK SYMBOL SYNC LOOP FOR A QASK INPUT 

In a QASK system, one has available two independent MASK signals as 
inputs to the symbol synchronizer. In particular, assuming perfect carrier 
synchronization, the output of the in-phase and quadrature carrier demodulators 
of the QASK receiver can be written, respectively, in the form 


14 


JPL Technical Memorandum 33-720 



yj(t) = Sjft, c) + nj(t) 
YqU) = s Q (t, e) + n Q (t) 


( 38 ) 


where 


Sj(t, «) 


a n p [t - nT - t] 


s Q (t, e) 


S b n P [t - nT - 0 


(39) 


Here p(t) is again the unit rectangular pulse defined by p(t) = 1, 0 < t < T, 
p(t) = 0; all other T, and the quadrature amplitudes a^ and independently take 
on equally likely values i& and j6 with i, j = ±1, ±3, • • • , ±(K-1). Also n^.(t) and 
n^(t) are uncorrelated (independent) white Gaussian noise processes with two- 
sided spectral density N q / 2. If the carrier synchronization is not perfect, 
then the carrier demodulator outputs are respectively given by [25, Eq, 7]. 

It is a simple matter then to modify the MASK symbol synchronizer of 
Fig. 1 to accommodate a QASK input. Fig. 16 illustrates such an implementa- 
tion. Basically, each input of (38) is fed through identical in-phase and quadra- 
ture circuitry with the resultant error signals added to produce the required 
symbol sync estimate. Thus, it appears at first glance that for the same signal 
step size 6, the QASK symbol sync loop has an S-curve whose values are twice 
that of the MASK loop. However, one must remember that the average trans- 
mitted power of a QASK signal set is S = [2(K 2 - l)/3]6 2 i. e. , twice that of an 
MASK set. Thus, for equal average transmitted powers, the p. d. f. of \ for the 

QASK - K^ symbol sync loop of Fig. 16 is identical to that of the K-level MASK 

2 

loop of Fig. 1. In the above, the notation QASK - K denotes a QASK system 
2 

with K (K X K) signals. Finally, one can relate the conditional error probability 
of a QASK - K^ system to that of a K-level MASK system by 


P E (X) 


QASK 


= 2P e (M 


MASK 


1 " 7 


2 p E<M 


MASK 




(40) 
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where P (X) is given by Eq. (34) with replaced by R^/2, Figs. 17, 18 

illustrate the average error probability performance of a QASK system obtained 

by substituting (38) and (40) into (37) and performing numerical integration on a 

digital computer. As in the MASK case, the average error probability is plotted 

2 

vs R, with b, as a parameter for £ = 1 and K = 4, 16 respectively. The case 
? d a 

K ~ 4 corresponds to quadriphase' signaling. 


G. CONCLUSIONS 

Summarizing, a symbol synchronizer for MASK and QASK signal sets has 
been proposed and analyzed insofar as its rms sync jitter performance and the 
effect of this jitter on system error probability performance. The sync jitter 
performance has been shown to be virtually independent of the number of signal 
levels, K. As in the binary case, the jitter performance improves as the 
window width £,T (the integration interval in the mid -phase arm) is reduced. 

The trade-off in reducing the size of the window is that the loop's cycle slip 
performance degrades. Thus, a compromise value of window width must be 
chosen. Typically, in the past (for binary systems), a quarter window (£ = l/4) 
synchronizer has been shown to yield adequate jitter and slip performance. The 
same conclusion would hold here in the multilevel case. From the results 
presented on error probability performance, we see that a symbol rate to loop 
bandwidth ratio 6^ = 250 (for K = 2 and £ = 1) and a 6^ - 100 (for K = 4 and £ = 1) 
are sufficient to reduce the signal-to-noise ratio degradation (relative to the 
ideal performance) to around .1 dB, In a practical system, 6^ would be on the 
order of thousands and £ would be less than unity; thus the practical performance 
of the symbol synchronizer would be virtually ideal. 

Although the phase detector characteristics (loop S-curves) obtained for 
tracking purposes have a rather jagged appearance, it must be remembered that 
for the system parameters of interest, the loop tracking range would be 
restricted to a small portion of the S-curve in the neighborhood of zero sync 
error where the S-curve can be assumed linear with a slope depending upon 
signal-to-noise ratio. Also, since these S-curves were obtained by averaging 
over all possible sequences of three symbols, each symbol ranging over K 
levels (the case of interest when tracking the input data), the S-curves for 
acquisition would be somewhat more regular in shape since there a specific 
acquisition sequence is transmitted. 


16 
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APPENDIX A 


EVALUATION OF THE LOOP NONLINEARITY g(k) 

To evaluate g(\) defined in the main text, and considering first the 
average with respect to the noise, we have' 1 ' 


n 


M = E n{-r[ 


3} 


Q k (c k + M k + N k ) - Q k (c k+1 + M k+1 + N k+1 ) 


( [M‘ + N 1 ] r- -} 

n | “2 L Q K ^ c k + M k + N k^ " Q K ( C k+ 1 + U k+ 1 + N k+ 1 J 


,{«< 


- E_ {4-\q Z k (c k + M k + N k ) - Q* (c k+1 + M k+1 + N k+1 )]j 


But (see Appendix E) 
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erf z = 
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e‘ t2 dt 


‘'Herein we omit the conditioning on X. in our notation, assuming it to be 
understood. 


JPL Technical Memorandum 33-720 


A-I 



Thus , 


the first term in 



is evaluated as 



where 



It follows from the independence properties of the noise random variables 
defined in the main text that 



Thus, the second term of reduces to 

I E n { M k °K [ c k + M k + N k] ' N i + 1 °K [ c k + 1 + M k+ 1 + N k + l] } 


In view of (A-l), this reduces to 
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where 


'(k- £/ 2)T + c 


Mi ' 
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n(t) dt 


' ( k - 1 / 2 ) T + e 


-(k-l/2)T+« 


N," = 
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n(t) dt (A-7) 


'(k-l + £/2)T + « 


Note: Mj^ 1 and N^ t are mutually independent and 


N 
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It can be shown that, 


|erf (a + bx)| = erf 


J' • ~ 2 2 


1 + 2b cr' 
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(A-9) 


Applying (A-7) - (A-9) to (A-6), we get 
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It can also be shown that. 
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Applying (A -11) to (A -10), finally gives „ 
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or in terms of <r , 


e r7fJ ex M^r exP] 


where (see Appendix E) 
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Finally, to evaluate the third term in ^ n | e j c ) P » we define {see Appendix E) 
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Thus, the third term of E n |ek|is evaluated as 
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It now remains to average je^j over the signal distribution. F 
the main text. 
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where d^ 4 a^/6 takes on values ±1, ±3, * • * , ±(K-1) with equal probability 1 /K. 
Substituting (A-18) and (A-19) in (A-17) and averaging overall possible equally 
likely sequences formed from the symbols d^_ ^ and d^. 
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where we have introduced the shorthand notation 


erf (j, m,£ ) 4 
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and 



(A - 22) 


After considerable algebraic manipulation {A- 20) may be simplified to 
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m j i 


1 < X <I 

2 K ~ 2 


where is the Neumann number, i, e, , € q = 1; = 2, H > 1, and the sums on 

m and j range over the values 1, 3, • *• , (K-l) while the sum on £ ranges over 
the values 0, 2, 4, • - • , (K-2), 
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APPENDIX B 


EVALUATION OF THE LOOP SPECTRUM, S(0, 0) 


We begin by considering 
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‘ We again omit the conditioning on X in our notation, assuming it to be 
understood. 
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where exp^ (-C /2(r^) and 0^ 7 (-C /2<r ) are given by (E-ll) and (E-18) respec- 
tively, then. 
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For m = I 
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TERM NO. 3 
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By comparison with Term No. 2, we write: 
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For m = 0 
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where from (E-27) and (E-28) 
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with cr' given in (14) of the main text, then 
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Term No. 10 = 


i 2 

it 


1 (4) 

K » v^< 


For m - 1 


+ 1 (4)/^±i\ - 2lf 2 >(4M l' 2> 


K 


\Jz cr 


K 


v/2 0-,7 K \Ji 


(B-35) 


Term No. 10 = 




1 (2)( C k \ 1 (2) [ C k+1 | _ 1 (2) / C k ] 1 (2 ) / C k+2 


K 


VzV,, / * 


K 


S*J K VA<r 
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TERM NO. 11 


£b 


E n | Q K (c k + M k + N k^ " Q K <C k+l + ^Si+l + N k+1* 


] 


[ c 


X Q„ (c, l + M. ^ + N ) - Q^. { Cl , + M , . + N. 

1 K ' k+m k+m k+m K k+m+1 . k+m+1 k+m+ 


4 
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For m = 0 


^ b k 

Term No. 1 1 = - 



Form = 1 


Term No. 11 = - 


TERM NO. 12 



For m - 0 


Term No. 12 = 









(B-39) 


C k +M k + N k>- Q K < C k +1 +M k+1 + N k+1>] 

+ M k+m + N k+m^ “ °K ^ C k+m + l + M k+m+l + N k+m+l ) ]| 

(B-40) 




(B-41 ) 
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For m = 1 


Term No. 12 = 


TERM -NO. 13 



For m - 0 


Term No. 13 = - 


For m ~ 1 


Term No. 13 = - 



^ C k + M k + N k^ “ Q K ^ C k+ 1 + M k+1 + N k+lU 

+ M k+m + N k+m^ " Q K ^ C k+m+l + M k+m+ 1 + N k+m+l 

(B-43) 





(B-44) 



exp K 
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TERM NO. 14 


8 E n j[_ Q K (C k+m + M k+m + N k+m } “ Q K (c k+m + l * M k+m + l + N k+m+l ) ] 


x [q k (c k + M k + N k ) - Q k ( c k+1 + M k+1 + N k +l^]| 


(B-46) 


For m = 0 


Same as Term No. 11 


(B-47) 


For m = 1 


Term No. 14 = - 


1^5 1<2) ( C k+l ) ,<2)(4k+2j ( C k ) (3) [ c k+l 

8 l K W K WJ K W K W,> 
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TERM NO. 15 


" ! E n j M k [°K (c k+m + M k+m + N k+m> " °K < C k+m+l + M k+m+l + N k+m+l>] 

* [°K (c k + M k + N k* " °K (c k+l + M k+1 + N k+ 1*] j (B-49) 


Same as Term No. 12 


(B-5G) 


Term No. 15 = - 


8 \Z2tt 


ex P K 


C k ) A2)[ c k +l ) (2) f C k+2 

^ K WJ K \yr-v 
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TERM NO, 16 


1 

8 ~n 


- 77 E <f N,' . , IQ?, (c. , + M, , + N. , ) - Q?, (c + M , , - + N. , , ,) 

K k+m k+m k+m K k+m+1 k+m+1 k+mtr 


| N k+l [' 

X K (C k + M k + N k> - °K < c k+l + M k+1 + N k+i>]} 


] 

(B-52) 


For m = 0 


Same as Term No. 13 


(B-53) 


For m = 1 


Term No. 16 



exp 


(3) 

K 



0 


<2M 

Kl 




(B-54) 


We must now average each of the terms (B-2) through (B-54), over the 
signal distribution. The procedure will be to collect all terms corresponding 
to a given value of m, and then average over all pos sible sequences where c^ f 
and are defined in (A- 18), (A- 19) and in addition 


^ = d k +m -2 X + d k + m-l f 1 


^ = d k +m -l X+d k +m <* 


for X < — 


(B-55) 


k+m _ , . 

6 k+m 
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'k+m 


= d. , 9 X + d. , , (1 - X) 

6 k+m -2 k+m- 1 


'k+m + 1 


d k+m- 1 ^ + d k+m (1 


- M 


£ 1 
for | < X < ^ 


(B-56) 


k+m 


d k+m - 1 ^ 


Summing first over the first nine terms, i.e. , (B-2) through (B-33), and evalu- 

vT^ 

ating the results at X = 0 gives 


For m = 0 




- 3 u ! + (K . u 2 

K - 1 



12*j 

2 2 

K (K - 1) 




The sum over m is for values m = 1, 3, 5, • • • , (K-l) while the sum over i 
is for values H = 0, ±2, ±4, • • • , ±{K-2), 
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Formal 


n, s { e k e k+l} _ 3 ^ R d 


E - -J_ e k e k+1 
0-4 

v ” d 


K 


? 2 

(K - 1) 


j[ y erf [ (m ■ 


m £ 


2 ? ?? 

KT (K - IT (K - 1) 

22 1 


K 2 (K 2 - 1) 


* JxE m E er£ [ <m - £, ^]E m 'E exp [■ (m ' - r > 2 R d 


m £ 


m ’ £ 


K Z (K 2 - 1) TT 


|E E exp [• (m • t)Z R d| 

y m £ * 




K (K - 1) 


E m2 E £ erf £(m - £ ) 


m £ 


X^r. E m E e exp [" (m ■ 1)2 R J 


K (K“ - 1) v ^ l 



ft (m - £ ) exp [- (m - £) 2 rJ 


Thus, for the first nine terms 


(B-58) 


E lef / + 2E -fc. c, 
n, s \ k; n, s j k k+1; 


ff 4 £R 


3 (K - l) 2 (K - l) 2 iR' 

7 ? • 

K 2 - I 2 


2 2 

K* (K £ - 1) 


R d E m E erf [ <m ■ e> 


EE-4 


(m - £) 


-■J 


m £ 




K {K - i) 


[ z + e R j ( !S ^r i )] EE'-'h- 1 ^] 
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(B-59) 

B-l 5 



Summing up the contributions of terms ten through sixteen, i. e. , (B-34) through 
(B-45) and evaluating at K = 0 we find that 


E 

n, s 



+ 2E |e e , / 
n, s l k k+lj 

u 4 ^Rj 

v d 


0 


(B-60) 


Combining (B-59) and (B-60) and recalling from (12) and (14) of the main text that 

o- 2 = N /2T, then 
v o 


h(0) A - 5( °j ° - 


K) K\ 

\2T/ \2T/ 


R(0, 0) + 2R(1, 0) 

2 ea d 


3 (K - l) 2 (K - l) 2 > , 

~ i?T7 1 s d 


K 2 {K 2 - 1)' 


R 


<z™z 


erf |(m - I } JR^ |+ 


ZZ exp j- (m - i ) 2 RU 


m £ 




K (K 


^ + erf [<m - I) JSj] 

* L ' ' -J O 


m £ 


where from ( 19 ) of the main text 


R(0,°) 4 E 



R(l, 0) 




(B-61) 
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APPENDIX C 


SLOPE OF THE LOOP S-CURVE AT THE ORIGIN 


The purpose of this appendix is to evaluate 


K 

g 


A d Sn (X) 
dX 


| X = 0 


(C-l) 


where g n (X) == g(X)/S with g(X) given by (11) of the main text. Differentiating 
(11) with respect to X and evaluating the result at X = 0 gives 


K 

g 





I 


" " 

(-j + m) exp 

-R' <-m + t ) 2 
a j 

+ (j + m) exp 

-R^(-m + l) 2 


j 


_ 



h -i 




r -l 

(-j - m) exp 

-R^(m + f) 2 


+ (2m - j) 

erf 

+ /) 








erf 

- 

^R^(-m + i) 


- (2m + j) 

erf 

■ 

J%(- m + J) 


- 
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C 



or 


K 

g 


2{K 2 



2exp 


r- 

-R^ (m + i ) 2 

L 


4- 2 exp 


- 




- 

-R' a {-m + f) 2 


+ 4m 

erf 

J R' a (m + f) 







1 


) 

J R' d (-m+ l) 



. 

J 

) 


(C-2) 


A-fter considerable algebraic manipulation, we arrive at the simplified, result 


K 


3K 


f 

\ 


E 


+ K 2 - U- erf (n Jh.) 
2 < K “-‘)I n=-K+3,“+5, K 


^ erf (nJ-R^) 

n=l,3, ' * • , K- 1 


E 


n 


=K+1, K+3, * * ' , 2K-3 


(ti - K 4- 1) 
K 2 


2 1 


erf 


(n /RU 
\ d 


i 

2N TT 


E 


n= -K+3 , -K+5, 


1 + 


n - 1 
K 


exp ( -n R’ ) 
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n= 1,3, * * * , K-Z 


exp { -n 2 R^) 


* E 

n=K+l, K+3, ' ' ' 2K-3 


1 


n - (K - 1) 
K 


exp ( -n 


2 



(C-3) 
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APPENDIX D 


DERIVATION OF THE CONDITIONAL ERROR PROBABILITY 
P E (X) FOR A MASK SYSTEM WITH A FIXED 
SYMBOL SYNC ERROR 


Consider the situation where the symbol sync is in error by 6 = \T sec. 
The cros s -correlator statistic on which a decision is made for the nth trans- 
mitted symbol now depends on the received data in the time interval 
(n - 1)T + t < t < nT + e where e is the sync estimate obtained from the 
symbol synchronizer. Stated mathematically, 


Cl 

n 

where is the nth detected symbol and y(t) = s(t, t) + n(t) is the correlation 
detector input with s(t f e) defined in (1) of the main text. Making the change of 
variables t 1 =t-(n-l)T-e in (D-l) and separating the signal and noise 
correlation components, we get 


nT+< 


= Q 


KIT 

( J(n-l)T+€ 


y(t) p[t - (n - 1) T - t ] dt 


(D-l) 


/\ 

a 

n 



Q 



(D-2) 


where 


Q s - Vl r s [- T <‘ - X >J +a n r s ^ T ) + a n+l r sf T < 1+X )] 
T 


D-3) 


Q = A 
n T 


p (t*) n[t' + (n - 1)T + £] dt’ 


and r (r) is the normalized signal pulse autocorrelation function defined by 
s 


T-t 


r s (T) = T 


P s (t') P s (t' - T) dt' = i 


P s (t' + r) p s (t') dt’ 


(D-4) 

D-l 


JPL Technical Memorandum 33-720 



Actually, if X is known to be positive (or negative), then only two symbols are 
involved, namely, the nth and one of its adjacent symbols. 

The average error probability conditioned on knowing X {assumed posi- 
tive) is given by' 1 ' 


p E m - 


E 

j, rn 


Pr ob \ 


error a 


n- 1 


= j5 , a = m6 


}Prob{a nl 


- j6, a 


n 


= m 6 


Since the input symbols occur with equal probability and are independent, 


P E (X) = Prob * m6|a n = m6, a nl = jo) 


J, m 


[Prob{Q n + j6R g [ -T( 1 - X)] 


E 

j m=±l, ±3, • ♦ • , ±{K-3) 

+ mSR g (XT) > (m + l)s| + Prob | Q R + j&R g (-T(l - X)] 


+ m6R g (XT) < (m - 1)6 j j 


+ A2^[Pr° b ) Q n + j6R s [-T(l - X)] + (K - l)6R s <VT) < (K - 2 ) 5 ) 


+ Prob | Q n + j6R s [ — T ( 1 - X)] - (K - l)6R g (XT) > - (K - 2)6 ) 1 


Since the random variable Q is Gaussian with zero mean and variance 

n 


(D-5) 


<r Q = E 


K} = 


N o s 

ZT 


(D-6) 


n 


Unless otherwise noted, the sums on m and j range over the values ±1, ±3, 
• • • , ±(K-1)« 
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then 


p e ( x) = zzY Y [? e H^[ (m + 1] - mr s (xT) 

j m=±l,±3,---,±(K-3) 

- jr s [ -T( 1 - M]]} + }erfc{ - 1) + mr s <\T) 

+ jr s [-T(l -M]]}] + -^J][lerfc{^[.(K-2) + (K- l)r s (VT) 

j 

+ jr g f-T( 1 - M]]}+ jerfc {m -(K - 2) + (K - l)r s (\T) 


- x «]} 


?E E 


erfc 




m + 1) - mr (\T) 


j 'm-±l, ±3, * * ‘ , ±{K-3) 


-jr s [-T(l -X)]]} + -^^erfc{/r k [-(K-2) 

j 

+ (K - 1 ) r g ( XT ) - jr s f-T(l - X)]]) 


(D-7) 


where R' is defined in (13) of the main text. Combining the two terms of (D-7) 
into one and considering now the answer for positive and negative we get' 1 ' 


p E (X) 


-iy er fc{/i^[< - (f - l)r s (|\|T) - jr s (T - |X|T)]} 
J > f 


(D-8) 


The summation on 


f ranges over the set of values 0, ±2, ±4, • * * , ±{K-2). 
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For rectangular pulses, i. e, , p(t) - 1, 0 < t <T, p(t) - 0; otherwise. 


1 - ^ 
rj» 


r (t) = 

s 


T 


< T 


| > T 


(D-9) 


Thus, for this case (D-8) reduces to 


P E (X > = - (I - 1) (1 - IM) - J IM]> 

K . / 

J. * 


(D- 10) 
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APPENDIX E 


EVALUATION OF THE MOMENTS OF VARIOUS FUNCTIONS 
ARISING FROM THE PASSAGE OF SIGNAL PLUS 
NOISE THROUGH A K-LEVEL 
QUANTIZER 

i. E n {Q K (c + y>} 

Consider the mean of the output of a K-level quantizer with transfer 
characteristic Q K (x) defined in (7) of the main text, when the input is signal 
plus noise, C + y . Letting y be a Gaussian r. v. with zero mean and standard 
deviation <r . then 

XT * 



(E-I) 

where the notation erf^-(x) is intended to represent a generalization of the 
standard error function erf x, defined in (18) of the main text to the K-level 
situation under investigation. When K = 2 and 6 = 1 , the generalized error 
function reduces to the standard error function, i. e. , erf^(x) = erf x. Substi- 
tuting (7) of the main text into (E-l) gives 



£ 

=±1, ±3, ' * ' , ±(K-3) 


■(H-l)6-C 


yfZTT (7 


exp I - 


y I)6-C 


2cr 


1 6 dy 


+ 




y .7(K-2)6-C 

-(K-2)6 -C 


exp j- | (K - I) 6 dy 

2o" 


\]Z TT C 


exp / - -L i r -(K - 1) 6 j dy 
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or 


erf 


K 


(T) 


E 


~ f 6 <; erf 


l=±l, ±3, • • ■ ,±(K-3) 


(i + 1) 5 - C 
\J~Z tr 


+ erf 


-(i - 1) 5 + C 
’ v^o- 


-A (K 


-1)6 erf 


(K - 2) 6 - C 

/2<r 


+ ^-(K - 1) b erf 


(K - 2) HC 


yfz tr 


(E-2) 


Finally, factoring out 6/2 gives 



U. E n {Q£<C + y)} 

We now consider the mean of the K -level quantizer output squared . If 
K = 2 and 6 = 1 t then the quantizer becomes a signum device whose output 
when squared has value unity regardless of the value of the input. Thus, the 
generalization of this result for the K-level quantizer is 


E-2 
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(E-4) 


E nK' c + y>) - 


f Zo y/ 


y ^-co 


exp |--^)Qk (C + y * dy 
2 cr 


/ 2 1 

The superscript notation on the function ( x ) refers to the fact that we are 

considering the mean of the square of Q t< -(x) . Also, from the above, for K = 2 

( 2 ) ^ 

and b = 1, we have that l^ 7 (x) ~ 1 • Squaring (7) of the main text and substi- 
tuting the result in (E-4) gives 


( 2 ) 

K 







y J (K-2}6 -C 
-(K-2)6-C 


exp | - — ~~\ [(K - 1)6] dy 
2<r 

y 


J*" (r - 


y J - as 


exp l-^y) [-{K - 1)6] 2 dy 
2o- 

y, 


(E-5) 


or by inspection. 
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It is a simple matter now to generalize the above to the case where the K-level 
qpantizer output is raised to an arbitrary power. Thus, 



(E-7) 


E-4 


JPL Technical Memorandum 33-720 



and erf^(x) 4 erf^(x) , Note that for K = 2, 6 = 1, and arbitrary m, 

K tv 


erf^^x) = erf x 


(E-8) 



Applying (A- 11) to (A-9) results in 
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where 


exp 


(tn) 

K 


and 



m 


6 '(K-l) m 


(C-{K-Z)6) 


^ 


E « 

£=±1 , ±3 , • * • , ±(K-3) 


m 


|exp j- 


(C-(£. 1)6) 

2<r 2 

y 


- exp 


[ (c - a + ns) 


2a- 


ill 


(E- 14) 



(E- 15) 


Also, for K = 2, 6=1, and arbitrary m, 



From (E-6) and the steps leading up to (E-10) we can immediately write 

> that 
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where 



For K = 2 and 6 = 1, 



More generally, for m even. 



(E- 19) 


(E- 20) 


where 





1) 


m 


- exp 


(C + (K - 2) 6 ) 
2 or 2 


+ exp 



^ [ eXP f (C ' ^ 

£=±1,±3,- .. ,±(K-3) L 1 y 
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we see that 
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exp 


( c - (J - 1 ) 6) 

2<r 2 

y 


+ -i- 



exp 




where erf-^(x) is defined in (E-3) and 
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where 



(x) is defined in {E-6) and 
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Fig. 2. Equivalent Data- Trans itio: 
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Fig. 4a. Normalized Loop Nonlinearity vs Normalized Symbol Sync Error; | - 0. 5 
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Fig. 14. Error Probability Performance of MASK System with Noisy 

Symbol Sync Reference; K=2 






Fig. 15. Error Probability Performance of MASK System with Noisy 

Symbol Sync Reference; K=4 







Fig. 16. Data- Transition Type of Symbol Synchronizer for QASK Signals 
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